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ABSTRACT 

The supermassive black hole at the Galactic center harbors a bound cluster of massive stars that should 
leave neutron-star remnants. Extrapolating from the available data, we estimate that ^1000 radio pulsars may 
presently orbit Sgr A* with periods of <100yr. Optimistically, 1-10 of the most luminous of these pulsars may 
be detectable with current telescopes in periodicity searches at frequencies near 10 GHz, where the effects of 
interstellar scattering are alleviated. Long-term timing observations of such a pulsar would clearly reveal its 
Keplerian motion, and possibly show the effects of relativistic gravity. We briefly discuss how pulsar timing 
can be used to study the dynamical and interstellar environment of the central black hole, and speculate on the 
prospects for astrometric observations of an orbiting pulsar. 
Subject headings: black hole physics — Galaxy: center — pulsars: general 
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1. INTRODUCTION 

Ten years of near-infrared observations of the Galactic cen- 
ter have revealed the proper motions of nearly two dozen stars 
within 0"5 of the compact radio source Sgr A*. Significant 
astrometric accelerations measured for 8 members of the Sgr 
A* stellar cluster point to a common center of gravity coin- 
cident with the posit ion of Sgr A*, and imply a central mas s 
of (3-4) x 1O 6 M (Ghez et all 12003 ri LSchodel et alJEool . 
Stars SO-2 and SO- 16 have the most compact orbits yet iden- 
tified, with respective periods of ~15yr and ~30yr, eccen- 
tricities of ~0.88 an d ~0.95, and comparable pericenter dis- 
tances of -100 AU dSchodel et alJl200l iGhez et ai1l2003bt 
lEisenhauer et all2003l) . If the central mass is confined within 
100 AU, the implied density is >10 16 M Q pc~ 3 , which essen- 
tially rules out existing models alternative to the hypothesis 
that the cent ral object is a supe rmassive black hole (BH; Maoz 
1998i lGhe7"et al.l2003htlSchodel et al 12001 . 

Evidence from the near-infrared spectrum of SO-2 
iGhez et alJ l2003a |) and th e inte grated spectrum within ~0."5 
of Sgr A* (benzel et al.lll997t jEckart. Ott. & Oenzeilll999t 
iFiger et all 120001 iGezari et all [20021) suggest that the ob- 
served Sgr A* stellar cluster is largely comprised of luminous 
(~10 4 L©), early-type (09 to B0) stars. If these stars are near 
the main sequence, their masses are 10-20M Q . How these 
stars came to reside so near the supermassive BH remains 
a puzz le; for discussions and references, see iGenzel et al.l 
(2003) and Ghez et al. (2003b). Nevertheless, the existence 
of a cluster of massive stars tightly bound to Sgr A* has im- 
portant implications. 

Stars of mass 10-20M Q have nuclear lifetimes of — 10 7 yr, 
and leave neutron-star (NS) remnants. Therefore, we expect 
a significant number of NSs to be bound to Sgr A* in orbits 
similar to those of the observed cluster stars, as well as in 
more compact orbits. Source confusion so far inhibits the dis- 
covery of stars with orbital periods of <10yr about Sgr A* 
( IGenzel et alJl2003tlGhez et al.l2003bl) . though we anticipate 
that massive stars and NSs populate this region. The most 
exciting possibility is if some of the NSs orbiting Sgr A* are 
detectable radio pulsars, an idea first considered in the pre- 
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scient article bv lPaczvnski & Trimble! dl979l) . In §0 we es- 
timate the total number of normal radio pulsars — i.e., those 
with surface magnetic field strengths of — 10 n -10 13 G — that 
may presently orbit the central BH with periods of <100yr. 

Radio-wave scattering in the interstellar plasma poses the 
largest obstacle to discovering pulsars near Sgr A*, where 
column the density of free electrons is very high. At observ- 
ing frequencies of ~ 1 GHz, pulsed emission from this vicinity 
suffers severe temporal broade ning, proh ibiting the detection 
of pulsars as periodic sources dCordes & Lazi J1991) . Rela- 
tively high frequencies of > 10 GHz are required to alleviate 
the effects of scattering. These issues are addressed in § |3j 
where we estimate the number of pulsars orbiting Sgr A* that 
may be detectable with current telescopes. 

The Keplerian motion of a pulsar orbiting Sgr A* would 
be clearly apparent in its long-term timing properties. Rela- 
tivistic gravity may introduce measurable deviations from the 
best-fit Keplerian timing solution, depending on the orbital 
parameters and timing precision. Various arrival-time delays 
and secular effects are quantified in § |4] In § |5j we consider 
what pulsar timing can teach us about the accretion flow onto 
the Galactic BH and the stellar dynamical environment, and 
investigate the possibility of astrometrically monitoring an or- 
biting pulsar. 

2. RADIO PULSARS ORBITING SGR A* 

The observed Sgr A* stellar cluster occupies the central 
~4000AU about the BH. Although the observational cen- 
sus is incomplete, this volume likely contains at least several 
dozen massive stars. However, it is not the present popula- 
tion of massive stars, but rather their predecessors, that are the 
progenitors of radio pulsars orbiting Sgr A* . Since the origin 
of the observed cluster stars is unknown, we can only spec- 
ulate on the history of the population of massive stars in this 
region. It is plausible that cluster stars are steadily or episodi- 
cally replenished as they evolve and leave NS remnants. This 
might be the case if the mechanism that feeds stars into the 
central ~4000AU is linked to the significant star-formation 
activity on larger scales. Observational evidence suggests that 
the central — 200pc of the Galaxy is a region of past and cur- 
rent star formation, with an average massive-star formation 
rate as large as ~10~ 3 yr , or ~10% of the rate in the entire 
Galaxy (Mezger et al. 1999; Launhardt, Zylka, & Mezger 
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2002). 

A plausible initial guess is that there are roughly as many 
pulsars within 0."5 of Sgr A* as there are massive stars. Thus, 
perhaps several tens of radio pulsars, with ages of <10 7 yr, 
orbit near the black hole. A simple extension of this argu- 
ment provides an estimate of the total number of active radio 
pulsars orbiting Sgr A*, including those much too faint to be 
detected by current telescopes. Assume that over a time of 
>10 8 yr an average number of ~4 0-100 NS progenitors or- 
bit Sgr A* with semimajor axes of <4000AU and periods 
°f Potb 100yr. A stellar lifetime of ~10 7 yr then implies 
a NS birthrate of ~10" 6 -10" 5 yr _1 . We further suppose that 
a large fraction of NSs turn on as radio pulsars shortly after 
birth. Radio emission terminates when a pulsar crosses the 
"death line" in the logP^-logfp plane, where P p is the pulse 
period (e.g., Rudak & Ritter 1994, and references therein). A 
coarse inspection of pulsar statistics in the logP ; ,-logf p plane 
shows that a median terminal age for normal radio pulsars is 
probably ~10 8 yr, depending on the mode l for the death line, 
and neglecting magnetic-field decay fe.g.. | Bh attacharva et al. 
H99l Itudak & Ritterlll994t iTauris & Konarll2001l) . There- 
fore, we predict that ~ 100-1000 radio pulsars presently orbit 
Sgr A* with P orb < lOOyr and ages of <10 8 yr. 

Typical NS "kick" speeds of v* ~ 100-300kms _1 are 
inferred from the proper motion s of ~100 isolated pul- 
sars in the Galactic disk fe.g.. [Hansen & Phinnevl 119971 
lArzoumanian. Chernoff. & Corde sll2002l) . Since the stars we 
are considering have orbital speeds of v or b > lOOOkms" 1 , an 
impulsive NS kick will usually cause only a small fractional 
change (~Vjfe/Vorb) in the orbital parameters. Therefore, the 
distributions of pulsar orbital parameters should be similar 
to those of their progenitors. For a stellar number density 
n(r) oc r -^ 1+ i)/ 2 about Sgr A*, the differential period distri- 
bution is p(Pcrb) oc Pjj if t he velocity di stribution is isotropic 
(e.g.,|Sch6del et al. 2003 ) . iGenzel et all J2003I) find that q ~ 
at <10" from Sgr A*, so that p(P OI b) is approximately flat. If 
there are ^1000 radio pulsars with P or b uniformly distributed 
over, e.g., 1-100 yr, then ~100 pulsars may have P or b < lOyr. 

3. PULSAR DETECTION 

Scattering of radio waves by fluctuations in the free- 
electron density causes angular broadening of radio images 
and temporal smearing of pulsed emission. Toward the Galac- 
tic center, observed scattering diameters of ~l"(z//GHz)" 2 
for Sgr A* and nearby OH masers imply a scattering timescale 
of ~300sO/GHz)"4 (lLazio & Cordeslll998albl) . where v is 
the observing frequency. Pulsars with P p < 1 s would thus 
be undetectable as periodic sources for v ~ 1 GHz. At v > 
5 GHz, the scattering time is <P p , but the flux density is re- 
duced according to the declining power-law spectrum, S v oc 
v~ a , followed by most pulsars. Cordes & Lazio (1997; see 
also Kramer et al. 2000) find that near Sgr A* the pulsed flux 
is maximized at v ~ 10 GHz for P p ~ 1 s, with a weak depen- 
dence on P p and a. At higher frequencies, the pulsed fraction 
of the flux is ~1, and the pulse duty cycle, e, approaches its 
intrinsic value; e ~ 0.05 is typical for normal pulsars. 

The minimum detectable flux density is 5 m ; n ~ 
CS„[e/(A^AWi„t)] 1/2 (e.g., IDewev et al.1 119851) . where 
C ~ 10 is the signal-to-noise threshold, S„ is the noise flux 
from the telescope and sky, Av is the bandwidth, N p = 2 is the 
number of polarizations, and f; nt is the integration time. For 
the 100-m Robert C. Byrd Green Bank Telescope 3 operating 



at v > 10 GHz, the noise flux is S„ ~ 25-50 Jy toward the 
Galactic center, and Av ~ 1 GHz. We then find that 

S mn * (20-40) M Jy(_y 1/2 (^L)- 1/2 . (1) 

From Green Bank, the Galactic center can be observed for up 
to ~8hr per day (S. Ransom, private communication), so that 
sensitivities of S m i n ~ 10/iJy may be possible. 

Observed Galactic pulsars have intrinsic 400-MHz lumi- 
nosities of (l-10 4 )mJyk pc 2 , with a cumulative distribu- 
tion f(>L4oo) ~L^L (e.g-lLyne. Manchester. & TavloJl985t 
Tavlor, Manchester, & Lynel ll993l iLvne et all 119981) . where 
L400 is in mJykpc 2 and is assumed to be <Cl0 4 . Pul- 
sar lum inositie s are typica lly defined by L„ = D 2 S V (e.g., 
iTavlor & ManchesteJI 19771) . where D is the distance. For 
a pulsar with S v oc v~ a at the 8-kpc distance of Sgr A* 
jEisenhauer et alJ 120031) . we have L400 — 64S„x Q mJykpc 2 , 
where x = z//0.4GHz, and S„ is in mJy. The distribution, 
p(a), of measured spectral slopes is roughly Gaussia n over 
a = 0-4.0, with a peak at a = 1.5-2 ( Lorimer et alJ fl995l 
iMaron et alJ2000)) . Utilizing the above pulsar statistics, along 
with the steady-state assumption of the last section, we now 
proceed to estimate the detectable fraction of pulsars near Sgr 
A*; ou r approach is quite similar to that of Cordes & Laziol 
(1997). 

At high frequenc ies, the d etection of shallo w-spectrum pul- 
sars is favored (e.g.. Lfohnston et alJl992tlWex. Gil. & Sendvkl 
1996). It is encouraging that ~10% of pulsars with a mea- 
sured spectrum have a = 0-1, although p(a < 1) is poorly 
constrained. We can crudely estimate the fraction, /(>5 m i n ), 
with flux densities greater than S m m by restricting to the range 
a = 0-1 and assuming a flat distribution, p(a < 1) = 0. 1 . After 
evaluating a simple integral, we find that 

«> s >=»(i^r (i^r < 2 » 

when i> 1. If p(a < 1) is taken to be a linear or quadratic 
function that vanishes at a = 0, each of which is consistent 
with the data, then /(>5 , m ; n ) — 2-3% at 10 GHz, and scales 
with frequency approximately as (ln;t:)" ( " +1) , where n = 1 (lin- 
ear) or n = 2 (quadr atic). If ^20% of pulsars a re beamed to- 
ward the Earth (e.g., Lvne & Manchester 1988), we conclude 
that ~1% of the active pulsars near Sgr A*, or as many as 
1-10 (see § |3, may be detectable with current technology. 
These pulsars would most likely have been born in the past 
several million years. The planned Square-Kilometer Array 4 
(SKA) will easily achieve sensitivities of <1 /iJy, and could 
perhaps ultimately find >100 orbiting pulsars. 

4. PULSAR TIMING AND RELATIVITY 

The dynamics of a pulsar orbiting Sgr A* are revealed 
through analysis of the pulse arrival times. The Newtonian 
dynamical signature of the supermassive BH should be clearly 
evident in a small segment of the orbit, via the acceleration of 
the pulsar. After a full orbit, a best-fit Keplerian timing model 
is obtained. Residuals of the Keplerian solution contain fur- 
ther dynamical information, including contributions from rel- 
ativistic gravity. In this section, we quantify various arrival- 
time delays and secular effects, and crudely assess their mea- 
surability. 
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TABLE 1 
Pulse Arrival-Time Delays 

Delay" Amplitude Width References 

Roemer b ~ 1 day Ml'Ipf' 3 sin i ~ 1 yrPj 1 

Einstein c ~1 hrM 2/ 5 3 P 1 1/3 e ~lyr^ 1 

First-Order Shapiro d ~30sM 6 . 5 1 ln[(l -e)(l -sin i')]| ~1 yrPi (1 -e) 3 / 2 (cos/) 1/2 1 

Second-Order Shapiro e -0.1 sA^ 6 ^P~ 2/3 (l -e)" 1 /cos/ ~1 yrP^l -efl 2 cos/ 2,3 

Frame Dragging f ~0.1 sMgj'p, 2 ^ 3 (1 -e) _ 'x/cos! — lyrP 1 (l-e) 3 / 2 ^cosi 2,3,4,5 



REFERENCES. — (1) Damour & Taylor 1992; (2) Dymnikova 1986; (3) Goicoechea et al. 1992; (4) Laguna 
& Wolszczan 1997; (5) Wex & Kopeikin 1999 

"The dimensionless variables used are M 6 . 5 = M BH /1O 6 5 M and P] = P orb /lyr. For simplicity, we have 
adopted w = 90° in estimating the amplitudes and widths. 

b Light travel time across the orbit. The Keplerian orbit is evident in the Roemer delay. 

c Combined effect of time dilation and the gravitational redshift. 

d Lowest-order relativistic propagation delay in the gravitational field of a point mass. 

c Next highest order contribution to the propagation delay that is independent of the BH spin. 

'Contribution to the net propagation delay due to the BH spin, in the special case where the spin direction is 
parallel to the orbital angular momentum of the pulsar. Here < \ < 1 is the dimensionless spin parameter, 
where x = 1 corresponds to an extreme Kerr BH. 



4.1. Time Delays and Secular Effects 

Each pulse arrival time, f, at the solar-system barycenter 
is related to the pulsar proper time, f', by f-fo = t' + Yli^i> 
where fn is a reference time, and the A,- are variable delays 
due to the Keplerian motion and relativity. The delays are 
functions of t' and the following Keplerian parameters: (1) 
BH mass, Mbh, (2) orbital period, P or b, (3) eccentricity, e, (4) 
inclination, i, and (5) longitude of pericenter, to. When frame- 
dragging is considered, we must also specify the magnitude 
and direction of the BH spin. Each delay may be characterized 
by an amplitude, the difference between the maximum and 
minimum values of A,, and a width, the timescale over which 
the delay shows the largest variation. 

In Table ^ we quantify the five delays with the largest 
expected amplitudes in the special case where u> = 90°, so 
that superior conjunction — when the pulsar is farthest behind 
the BH — coincides with pericenter passage. This simplifies 
the analysis and illustrates roughly the dependence on eccen- 
tricity. Dimensionless variables used in Table 1 are Me.5 = 
M BH /1O 6 - 5 M and P { = P Qlb /l yr. Not included in Table 1 are 
the d elays due to aberration o f the beamed pulsar emission 
(e.g., ISmarr &~B landford 19761) or the bending of light rays 
in the gravitational field of th e BH (Doroshenko & Kopeikin 
1 19951 IWex & Kopeikinlll999l) . Each of these delays has an 
amplitude of <1 ms for a wide range of P or b, e, and i. 

Several secular processes cause changes in the orbital 
elements over long timescales; these include the follow- 
ing. Emission of gravitational radiation by an orbiting pul- 
sar causes the semimajor axis to shrink on a timescale of 

~10 I3 yrM~^Vf /3 (l-e 2 ) 7 / 2 re.g.. lTavlor & Weisbergll 989h . 
which is typically too long to be of interest. The geodetic pre- 
cession rate o f the pulsar spin axis is ~0.° 1 M 2 J^P, 2/ ^(l-e 2 ) -1 
per orbit (e.g., Weisberg. Rom ani. & Tavlorfl9 89). This is ev- 
ident as a long-term change in the pulse profile, and its mea- 
surement depends on, among other things, the precise geome- 
try of the pulsar beam. For e > 0.9 and Pi ~ 1 it is conceivable 
that geodetic precession could be detected after several orbits. 
We now estimate the Newtonian and relativistic contributions 



to the secular apsidal precession of the pulsar orbit. 

Suppose that surrounding the BH is a spherically symmetric 
distribution of matter, in the form of mostly low-mass stars 
and compact objects (e.g., § 15. 2> . For a density profile p oc 
r~ 7 , the extended mass enclosed within a radius r is M e (r) oc 
r 3 " 7 (7 < 3). If M e (a)/MBH <C 1 for an orbit with semimajor 
axis a, the Newtonian contribution to the change in uo per orbit 
is 5 



Acjn : 



M e (a)(l-e 2 ) 3 - 7 



BH 



2tt 



(l+eCOS(/)) 3 - 7 



(3) 



where <f> is the true anomaly for the unperturbed elliptical tra- 
jectory. If 7 = 2 (a plausible choice) the above integral is ana- 
lytic, and we have 



AlJn = — 27T 



M e (a) 1-e 2 
M RH e 2 



1 



(l_ e 2)l/2 



1 



(4) 



where the minus sign indicates retrograde precession, which 
is generally the case for 7 < 3. For example, if 7 = 2, 
M i ,(a)/M BH = 0.01, and e = 0.9, we find that Aw N ~ 1°. 

The net apsidal precession rate also includes two relativistic 
contributions. In the Schwarzschild spacetime, the prograde 
advance per orbit is (e. g. JWeinbergl 1 97% 



Au> s 



-2/3, 



(5) 



For a spinning BH, f rame dragging introduces an addi- 
tional contribution (e.g., Jaroszvns kill998bHWex & Ko peikin 
1999): 

Aw FD =s -27"M 6 . 5 P;\1 - e 2 r 3 / 2 X cos^ , (6) 

where ip is the angle between the angular momentum vectors 
of the BH and the orbit, and < x < 1 is the dimensionless 
BH spin parameter. 

5 We determined the precession rate by computing t he orbit-av eraged 
rate of change of the Laplace-Runge-Lenz vector (e.s.. IGoldsteinlll980D . 
e = (GMBHr'v X h-r/r, where h = r X v. Weinbera 1 1972) uses the same 
technique to calculate the general relativistic precession rate. 
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4.2. Remarks on Measurability 

The degree to which different contributions to the tim- 
ing residuals can be resolved depends on the precision and 
number of measured arrival times. Typical precisions are 
St ~ (10" 3 -10" 2 )P p , or -l-10ms for P p ~ Is. If one aver- 
age arrival time is measured each day the pulsar is observed, 
and N ~ 100 arrival times are measured per orbit (P 01 -b > 1 yr), 
then we expect a root-mean-square (RMS) timing precision 
of e ~ St I y/N < 1ms over one orbital period. However, vari- 
ous potential sources of error may limit the net timing preci- 
sion in a single orbit to ~P P . In particular, stochastic "timing 
noise," which is largest for the youngest, most luminous pul- 
sars, may introduce net re siduals of >0. \P p after several years 
(e.g.,|Arzoumanian et al. 1994). In general, a precise assess- 
ment of the measurability of arrival-time delays and secular 
effects requires detailed simulations that cover a large param- 
eter space. Here we present some simple statements regarding 
the detection of the delays listed in Table 1, as well as apsidal 
precession. 

The amplitudes of the Einstein and first-order Shapiro de- 
lays given in Table 1 suggest that these effects should be easily 
measurable over a wide range in orbital parameters. This may 
indeed be the case for the Shapiro delay, allowing for an inde- 
pendent determination of sin/, a s lon g as perturbations to the 
orbit by stellar encounters (see § I5.2> do not have a significant 
impact. However, the Einstein delay can only be measured 
if the orbit undergoes ap preciabl e apsidal precessio n (e.g., 
iBlandford & Teukolskvll976tlDamour & Tavlorll99l . In ef- 
fect, extraction of the Einstein delay requires that the change 
in cos a; be resolved sufficiently, which is generally more diffi- 
cult than measuring ui alone (see below). For a characteristic 
RMS timing precision of e < P p , the second-order Shapiro 
and frame dragging delays may be just at the threshold of de- 
tectability for Pi ~ 1, unless the orbit is highly inclined and 
eccentric. 

The rough scali ng of the fractional error in th e measured 
value of uj is (e.g., Blandford & Teukolskv 1976) 



5d> 



10" 



kc 



NarbeAud sin/ 



(7) 



where £ is a dimensionless factor that depends on the initial 
values of the Keplerian parameters, N OI b is the number of or- 
bits over which the pulsar is monitored, Aud is the net apsidal 
advance per orbit in degrees, and e is in seconds. Even un- 
der rather unfavorable circumstances, where, e.g., k ~ 100, it 
might be possible to obtain a < 1 0% measurement of u> in only 
two or three orbits if Auod ~ 1 , thus facilitating the detection 
of the Einstein delay. However, the various contributions to u> 
cannot be determined independently. The degeneracy could 
be broken if two or more pulsars are monitored over several 
orbits; such an ambitious project probably must wait for the 
SKA. 

5. DISCUSSION 

Here we address three additional topics pertaining to obser- 
vations of a radio pulsar orbiting Sgr A*. We first discuss how 
pulsar timing can be used to probe the physics of the accretion 
flow onto the BH. This is followed by a short investigation of 
the effects of gravitational interactions between an orbiting 
pulsar and the surrounding cluster of stars and remnants. Fi- 
nally, we consider the prospects of astrometrically monitoring 
an orbiting pulsar. 



5.1. The Interstellar Plasma Around Sgr A* 

Timing observations of a pulsar orbiting Sgr A* can be used 
to derive the properties of the local interstellar plasma. As 
the pulsar moves, photons trace many different lines of sight 
through the plasma to the observer. A gradient in the free- 
electron density on the scale of the pulsar orbit introduces dis- 
persive and refractive pulse arrival-time delays that vary over 
the orbital period 6 . Each of these effects has a characteristic 
frequency dependence, pointing to a need for multi-frequency 
observations. Because the detection of shallow-spectrum pul- 
sars is favored at high frequencies (see § [3}, there is a good 
chance that a pulsar detected at 10 GHz will also be detected 
at 15-20 GHz for the same sensitivity. 

Chandra observ ations of the Galactic center 
(Baganoff et al. 2003) indicate that the density and tem- 
perature of the electrons at ~1" (~8000AU) from Sgr A* 
are n e ~ 100cm" 3 and kT e ~ 1-2 keV. The gravitational 
potential of the B H exc e eds th e thermal energy of the 
plasma inside the Bondi (1952) radius, Rb ~ GMbu/c 2 , 
where c s ~ {kT e /m p ) l l 2 is the thermal speed, assuming 
equipartition between electrons and protons, where m p is the 
proton mass. For the plasma around Sgr A*, we find that 
Rb ~ 10 4 AU. Within Rb, n e (r) depends on the physics of 
the accretion flow. We adopt n e (f) ~ 1 2 cm" 3 ( 1 4 AU /r)P, 
where models predict ~ 1—1.5 (e.g., iMelia & Falckdl2001l: 
I Yuan. Ouataert. & Naravanl2003l) . 

The plasma frequency is v p = (n e e 2 / nm,,) 1 I 2 ~ 
90kHz(10 4 AU/V)^/ 2 for the density profile given above. The 
index of refraction of the plasma is £(i/) ~ 1 - (v p /v) 2 /2, 
for v 3> v p . Radiation at different frequencies propagates 
through the medium with different group velocities of c^{v), 
so that the arrival time of a pulse is frequency dependent. 
The difference in arrival times at frequencies v\ and vi is 
proportional to {v\ 2 - v^DM, where DM is the dispersion 
measure, the column density of free electrons along the path 
of the pulse. Modulation of DM over the orbital period of 
a pulsar bound to Sgr A* can thus be used to constrain the 
density profile of the plasma within the orbit. 

Refraction due to the large-scale gradient of the electron 
density causes a net angular deflection of individual pulses 
that reach the observer. Consequently, there will be a vari- 
able geometrical time delay as the pulsar orbits. We assume 
that a light ray is refracted impulsively as it passes its clos- 
est approach, b, to Sgr A* — the "scattering-screen" or "thin- 
lens" approximation. The deflection angle is 9 P ~ B[v p {r = 
b)/i/] 2 /2, pointing away from Sgr A*. From the thin-lens 
geometry, the resulting excess propagation time, compared 
to that of a straight-line path, is ^r8 2 /2c oc v~ A , where r 
is the orbital radius. For an orbital period of 10 yr, incli- 
nation of / = 80°, and eccentricity of e = 0.9, the minimum 
possible impact parameter is a(l-e)cos/ ~ 10 AU, where 
fp(10AU) ~ 16MHz for (3 = 3/2. At v = 10 GHz, the ampli- 
tude and width of the refractive delay are, respectively, ~0.1 s 
and ~10 days. The strong frequency dependence of the re- 
fractive delay distinguishes it from dynamical effects. 

5.2. Dynamics of the Sgr A* Cluster 

A given orbiting pulsar will interact gravitationally with 
an unknown number of normal stars, other NSs, white 

6 Here we are considering the gradient in the spatially averaged electron 
density near Sgr A*. This is distinct from the small-scale turbulent density 
fluctuations that are responsible for angular and pulse broadening (see §|5}. 
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dwarfs, and stellar-mass BHs. An important prediction of 
our work is that thousands of NSs born over the past ~1 Gyr 
may orbit Sgr A* with periods of <100yr. In addition, 
Miralda-Escud e & Gould! (120001) suggest that ~10 4 stellar- 
mass (~10M©) BHs may have migrated due to dynamical 
friction to the central ~1 pc about Sgr A*; extending their re- 
sults, we find that perhap s ~10 2 suc h BHs may reside in the 
central 1 ".|Ch aname & Gould 12002) have discussed the pos- 
sibility of using positional information for ^50 millisecond 
radio pulsars within a few parsecs of Sgr A* to indirectly in- 
fer the presence of a large population of stellar-mass BHs. 

For a star orbiting Sgr A*, many weak gravitational pertur- 
bations accumulate over a relaxation time to yield a signif- 
icant net change in the orbital parameters 7 . In the absence 
of resonant effects, it is straightforward to show that near 
Sgr A* the relaxation time is r rel ~ (10/V S ) -1 (M B}1 /M s ) 2 P olb , 
where M s and N s are the the typical mas s and total nu mber, 
respe ctively, of the perturbing stars (e.g.. lRauch & Tremaind 
1996). If M B n/M s ~ 10 6 and N s = 10 3 -10 4 , we find that 
T r d ~ (10 7 -10 8 )F n rh- Resonant angular- momentum relaxation 



(Rauch & Tremaine 1996), which causes variation in only the 
orbital eccentricity and orientation angles, can act on a much 
shorter timescale of r re i, res ~ (M Bli /M s )P orh ~ 10 6 P orb . 

Random-walk fluctuations in the orbital energy and angu- 
lar momentum lead to a root-mean-square fractional change 
in one orbit of approximately (P ib/ r ) 1 ' /2 ~ 10~ 3 — 10~ 4 , where 
t is either the conventional or resonant relaxation time. The 
resulting change in a/c may be ^100-1000 light seconds for 
^orb ~ lOyr. Gravitational perturbations by stars and remnants 
may have measurable consequences for the timing analysis of 
a pulsar orbiting Sgr A*, possibly even masking the impor- 
tant relativistic time delays discussed in the last section. With 
the SKA it may be possible to obtain a large sample of reg- 
ularly timed, orbiting radio pulsars. If evidence of random 
gravitational encounters is found in the timing properties of a 
significant fraction of these pulsars, important clues could be 
extracted regarding the number, masses, and velocities of the 
Sgr A* cluster members. 

5.3. Radio Astrometry 

It would be an extreme challenge to image a pulsar or- 
biting Sgr A* with an interestingly short orbital period. 
The projected apocenter separation is less than a(\+e)/D = 

1/3 2/3 

20masM 65 f > 1 (1+e). Very long baseline interferometry 
(VLBI) would be required to resolve a faint pulsa r next to the 
very bright Sgr A* (-1 Jy from 1 to 10 GHz: lMe"lia & Falckd 
1200 ll) . Such observations would have to be conducted at 
v > 10 GHz to sufficiently reduce the scattering diameter 
of Sgr A* (see § |3}. With a current VLBI sensitivity of 
~lmJy, probably no orbiting pulsars would be detectable 
(see § |3- The SKA will have the requisite sensitivity, and 



the current design specifications call for an angular resolu- 
tion of ~ 10 mas at 10 GHz, and, of course, higher astrometric 
resolution. Astrometric precisions smaller than ~0. 1 mas at 
^ 10 GHz may not be attai nable from the surface of the Earth 
(e.g.,|(3iatteriee et al] l2004l) . The advent of highly sensitive 
space- or Moon-based VLBI instruments is then a necessary 
step toward <10/ias astrometry for Sgr A* pulsars. 

If Aujd is the net apsidal advance per orbit in degrees, then 
the angular shift of the apocenter position in one orbit is 



A9 a < 0.34masAw d M 6 . 5 f > 1 2/3 (l+e) 



(8) 



For a large eccentricity and Pi ~ 10, A9 a could be ~1 mas. 
It is at least conceivable that apsidal precession could be de- 
tected from the ground after several orbits. 

If the orbit is inclined with respect to the BH spin, frame 
dragging causes the line of nodes, fi, on the plane of the sky 
to advance in one orbit by an amount (e.g., Jaroszvnski 1998b) 



Ail 



FD 



i9"M b . 5 P[ [ (l-e 2 T 3/2 X smiP 



(9) 



7 This is due to the granulari ty of the stellar distribution. Long-term orbital 
} 14.11 is due to the smoothed potential of the 



precession, as discussed in 



where xsinV 1 ls the projection of the spin onto the orbital 
plane. The corresponding apocenter shift is expected to be no 
larger than a few tens of /ias. If the contributions from frame 
dragging to apsidal (see § and nodal precession are mea- 
sured, it would be possible to determine the magnitude and 
direction of the BH spin, but this is a highly unlikely prospect. 

The effects of gravitational lensing by the BH may be im- 
portant when the pulsar is near superior conjunction. For lens- 
ing by a point mass, two images are produced, one inside 
and one outside the Einstein radius, 9 E ~ {AGM^r / c 2 D 2 ) x I 2 , 
where r -C D is the orbital radius near superior conjunction. 
For r = 100-1000 AU, we find E ~ 0.5-1. 5 mas. Since 9 E is 
small, we consider here only the image outside 9 E . At superior 
conjunction, the angular separation between the pulsar and its 
imag e is 89 ~ 9 2 ,D/b 1 wh ere b = r cos i is the impact parameter 
(e.g., Jaroszvnski 1998a). We then find that 59^20 /ias/ cos/, 
which is ~1 mas for i = 89°. Therefore, gravitational lensing 
has a small astrometric signature in typical situations. 
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